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Abstract. Recent work on purely transmitting 'jump-defects' in the sine-Gordon 
model and other relativistic field theories is extended to non-relativistic models. In 
all the cases investigated the defect conditions are provided by 'frozen' Backlund 
Sh ' transformations and it is also shown via a Lax pair argument how integrability will be 

preserved in the presence of this type of defect. Explicit examples of the scattering of 
solitons by defects are given, and bound states associated with 'jump- defects' in the 
nonlinear Schrodinger model are described. Although the nonlinear Schrodinger model 
provides the principal example, some results are also presented for the Korteweg de 
00 \ Vries and modified Korteweg de Vries equations. 
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. 1. Introduction 

The study of impurities, or defects, has a lengthy history especially in the context of 
condensed matter physics (see for example the lecture notes by Saleur [1]). In the 
context of integrable field theories of various types there has been interest in studying 
defects both from a classical and a quantum point of view, either theoretically or from 
the applications perspective. In the quantum domain there are issues surrounding the 
extent to which integrability is compatible with the reflection and transmission typical 
of a defect. The pioneering work of Delfino, Mussardo and Simonetti [2] represents 
one point of view and more recent work of Mintchev, Ragoucy and Sorba [3] provides 
another. Typically, adding a 5-impurity to a classically integrable nonlinear field theory 
destroys its integrability, although there may nonetheless be interesting phenomena 
associated with defects of this type (see [4] for the behaviour of solitons in the sine- 
Gordon model with a 5-impurity, or [5] and [6] for studies of some aspects of the 
nonlinear Schrodinger model with a 5-impurity) . On the other hand, there may be 
circumstances where a different type of defect is able to preserve the property of classical 
integrability and it is an interesting question to investigate what those circumstances 
might be. Some years ago, it was noticed that several types of relativistic integrable 
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field theory permit discontinuities (a type of defect) without the property of classical 
integrability being destroyed [7]. Among these are free fields, Liouville theory, the 
sine/sinh- Gordon model and a variety of affine Toda field models (though possibly 
not all of them [8]). Besides preserving integrability, it was found that the conditions 
defining the defect allowed not only energy conservation (as would be expected) but also 
the conservation of a generalised momentum, including a contribution from the defect 
itself (which was unexpected because of the evident loss of translation invariance). An 
integrable discontinuity of this type will be referred to as a 'jump-defect' in order to 
distinguish it from a (generally non-integrable) 5-impurity 

The integrable jump-defects have a Lagrangian description and their integrability is 
ensured by the existence of suitably adapted Lax pairs. The specific form of the jump- 
defect conditions is quite striking and makes novel use of Backlund transformations 
which are well-known to be a bulk feature of each of the models considered (for example, 
see [9]). With a single jump-defect the setup is is quite straightforward and easily 
described. One of the simplest examples is provided by the sine-Gordon model, as 
follows. 

The sine-Gordon model in the bulk [10, 11, 12] is specified by the Lagrangian density 

1 171 ^ 

C=-(u 2 t -ul)-—{l-cos(3u). (1.1) 

A single jump-defect placed at x = is then described by altering the Lagrangian in the 
following manner, denoting the field on the left of the defect by u and the field on the 
right of it by v. The full Lagrangian consists of pieces from the bulk regions (x < and 
x > 0), together with a delta function contribution at x — 0. Thus, the new Lagrangian 
density is given by 

£ = 9{-x)C u + 0(x)C v - S(x) f ' 



(uv t - vu t ) + B(u,v) 



(1.2) 



with 



„ 2ma n fu + v\ 2m (u - v\ 

The form of the additional piece is required by integrability and leads to defect conditions 
linking the two fields u and v and their derivatives at x — 0, the position of the defect. 
The usual variational principle reveals that the defect conditions constitute a Backlund 
transformation 'frozen' at the defect. There are many interesting features of these 
defects, both from a classical field theory point of view - especially with regard to the 
behaviour of solitons - and within the quantum context [13]. There is no obstacle to 
having several defects at different locations (with the same or different parameters). 
Moreover, they can move and scatter amongst themselves [13]. As mentioned already, a 
novel feature exhibited by these defects is the manner in which they may exchange both 
energy and momentum with the fields on either side of the defect location. There is 
nothing surprising about the energy since time translation is unbroken, but the existence 
of a conserved - though modified - momentum is quite surprising since translation 
invariance is certainly lost. In fact, the defect potential (1.3) may be regarded as being 
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determined by demanding that it be possible to find a conserved generalised momentum 
functional. The ability to exchange momentum and energy with the fields on either side 
of it is a defining feature of a jump-defect, and will be used as a tool later in this paper. 

As a further example, and to make clear the distinction between a jump-defect and 
a 5-impurity, consider two free-field situations. The jump-defect can be described as a 
small field approximation to (1.2, 1.3). In other words, 

vf (1.4) 

= \{v 2 t -v 2 x -m 2 v 2 ). (1.5) 

On the other hand, in the same notation, a ^-impurity is described by the quadratic 
Lagrangian, 

C s = 9(-x)C u + 9(x)C v -8(x)- [auv - (u x + v x )(u - v)\ , (1.6) 

leading to free Klein-Gordon equations in the bulk together with the defect conditions 

u = v , v x — u x = au, x = 0. (1.7) 

It is not difficult to check by explicit calculation that the former permits a conserved 
total momentum functional, including a defect contribution, while the latter does not. 
There is no requirement for the linear jump-defect to satisfy u = v at x = 0, and in 
general there will be a discontinuity. It is easy to check that the linear jump-defect is 
purely transmitting, while the 5-impurity both transmits and reflects. The jump-defect 
does not possess a classical bound state but the 5-impurity does for a suitable range of 
a. 

Besides the sine/sinh-Gordon model there are other integrable equations which have 
arisen naturally in special physical systems. Most of these are non-relativistic but might 
nevertheless allow discontinuities related to jump-defect conditions of frozen-Backlund 
type. One such example is the nonlinear (cubic) Schrodinger equation (NLS); others 
are the Korteweg-de Vries (KdV) and modified KdV (mKdV) equations. 

As already mentioned briefly, defects in the context of integrable field theories 
have been discussed before in the quantum domain, starting with Delfino, Simonetti 
and Mussardo [2], and elaborated subsequently by others, including Konik and LeClair 
[14], and Castro- Alvaredo, Fring and Gohmann [15]. As a consequence of their work 
it appeared that a defect could not allow simultaneous reflection and transmission 
while maintaining integrability (in the sense that the reflection and transmission 
matrices together with the bulk S-matrix should satisfy a set of algebraic compatibility 
requirements) unless the bulk scattering matrix was independent of rapidity. On the 
other hand, the sine-Gordon jump-defect is purely transmitting as far as its behaviour 
with respect to solitons is concerned and evidence was gathered in [13] to support the 
idea that the pure transmission matrix discovered by Konik and LeClair really describes 
the quantum version of the sine-Gordon jump-defect. Recently, Mintchev, Ragoucy and 
Sorba have suggested alternative compatibility relations that might be satisfied between 
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the transmission and reflection factors and the bulk S-matrix [3]. Using this framework 
5-type impurities have been explored in the context of NLS [16]. 

The purpose of this paper is to demonstrate that jump-defects can be incorporated 
into the NLS model very naturally without spoiling integrability. Following an analysis 
of a few of the conserved quantities - actually sufficient to determine the form of 
the defect conditions - an argument based on a generalised Lax pair is given in 
sufficient detail to demonstrate how integrability is preserved. Similar but less complete 
arguments are also given in the context of KdV and mKdV. In all cases, remarks are 
made concerning the behaviour of solitons as they encounter a jump-defect. To date, 
the only genuine surprise relative to the sine-Gordon case is provided by the strange 
(and still mysterious) behaviour of a 'fast' soliton in KdV. It remains to be seen what 
the quantum version of the NLS jump-defect will turn out to be and how it will relate 
to earlier work, if at all. 



2. The nonlinear Schrodinger equation 

The nonlinear Schrodinger equation with a cubic interaction term will be taken to be 
defined by the field equation [11] 

iu t + u xx + 2u{uu) = . (2.1) 

This may be derived in the bulk using an action principle based on the Lagrangian 

C = - (uu t — u t u) — \u x \ 2 + (T 2 |w| 4 . (2.2) 

A (real, positive) coupling constant a 2 has been added but then scaled away in the field 
equation by redefining u. Henceforth it will be ignored. The sign of the cubic term is 
important for some considerations. For example, the sign chosen in (2.1) is appropriate 
for a model possessing soliton solutions (corresponding to an 'attractive' interaction). 

If there is a defect at x = then the bulk fields to either side of it will be denoted 
u and v, and a boundary contribution B, presumably depending on u,v,Ut and Vt (and 
possibly spatial derivatives), will need to be added. In other words, the full action will 
be 

r 



//*u /*oo 
dt / dx£(u)+B+/ dxC(v) 
J-oo JO 



(2.3) 



The corresponding defect conditions at x = are: 

OB_d_dB (2 A) 

du dt du t ' dv dt dv t ' 

with similar expressions for the conjugate fields. Note that (2.4) would need to 
be modified if the defect part of the action depended additionally upon the spatial 
derivatives; for the time being it will be assumed these are absent. 
The analogue of energy for NLS is the density 

£=\u x \ 2 -\u\\ (2.5) 
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and, since the system remains time-translation invariant despite adding the defect, 
the total energy including a defect contribution will be conserved. In the bulk, as a 
consequence of space-translation invariance, the momentum has the density 

V = i (uu x - u x u) (2.6) 

and it is certainly conserved. However, when there is an impurity the momentum is 
not expected to be conserved since translation invariance is broken. However, as has 
been noted before in other cases, for example in the sine-Gordon model, this is not 
necessarily so. Adding a jump-defect does not spoil momentum conservation since this 
kind of defect can exchange both energy and momentum with the fields u and v to either 
side of it. To demonstrate this, note that the total contribution from the fields u and v 
to the momentum density is given by 

/0 roo 
dx i (uu x — u x u) + dxi (vv x — v x v) (2.7) 
-oo JO 

and it easy to see that the time derivative of P will depend critically on the defect 
conditions since 

P t = (2\u\ 4 + 2u x u x + i(uu t - utufj - (2\v\ 4 + 2v x v x + i(vv t - v t vfj . (2.8) 

To allow conservation of this charge requires that the right hand side of (2.8) should (if 
possible) be a total time derivative of a functional of the fields evaluated at the jump- 
defect. Put alternatively, if the momentum, suitably modified, is to be preserved then 
the defect conditions ought to be chosen to ensure it. This is not quite straightforward 
to achieve (and the conditions associated with a 5-impurity do not have this property) 
but, bearing in mind (2.4) the following suggestion works perfectly. It will be seen later 
that it also fits more generally with the idea of integrability. It is sufficient to take 



B = fl 



2 dt 



u — V 
u — V 



where 



B 



1/2 | 

-[a -\u- 



+ B 



u\ + \v\ 



Q=(a 2 



\u 



2 ) 



1/2 



(2.9) 



(2.10) 



and a is a real parameter. Although this is not a completely transparent choice, it is 
based on knowledge of the NLS Backlund transformation and experience with the sine- 
Gordon equation. It will become clear that (2.9) has nice properties. With the choice 
(2.9, 2.10) the defect conditions (2.4) at x = become: 



u. 



i (ut ~ vt) 

n 

i (ut - vt) 



+ (u + v)tt + 



[u — v)(\u\ 2 + \v\ 2 ) 

n 

[u — v ){\u\ 2 + \v\ 2 ) 



and the time derivative of the momentum simplifies to: 

P t = ^ [UV ~ VU) X=0 , 



(2.11) 



(2.12) 
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and thus the quite attractive combination 

P — \{uv — vu) x=0 (2-13) 

is conserved. The requirement that the time-derivative of the momentum should turn 
out to be expressible as a functional of the fields evaluated at the defect is actually very 
strong, and finding an expression which satisfies (2.8) severely limits the choice of defect 
condition. 

On the other hand, as mentioned above, the total energy should be conserved 
whatever the defect condition is, provided it does not violate time translation invariance. 
Indeed, the total energy satisfies 

E t = ^(nB) x=0 , (2.14) 

implying that E — (QB) x= o is conserved. Notice that the discontinuity [v — u] x=0 at the 
jump-defect should not be too severe. In fact 

\u-v\ 2 x=0 < a 2 , (2.15) 

otherwise the energy would fail to be real and simply lose its meaning. In the limit 
a — > 0, the discontinuity disappears. In this sense, the parameter a controls the height 
of the jump. 

There is also a 'probability', or 'number', density M = uu which satisfies 

■A/i = i (u x u — uu x ) x . (2.16) 

Before checking that the NLS equation remains integrable even after adding a jump- 
defect it is also worth examining this charge in some detail. On the whole line this 
'number' is certainly conserved (assuming suitably decaying u at x — ±oo, or periodic 
boundary conditions) and it is a consequence (via Noether's theorem) of the continuous 
U(l) symmetry of the action under the constant change of phase 

u^e iA u A x = = A t . (2.17) 

When a jump-defect is added, it is natural to define, 

dxAf(u) + / dxAf(v), (2.18) 

-oo JO 

the defect condition (2.11) leads to 

N t = i [u x u - uu x ) x=Q - i {v x v - vv x ) x=0 = — 
It is clear from this that the combination 

1/2 



(2.19) 



x=0 



N + tt = N(u) + N(v) - (a 2 - \u - v\ 2 ) 



x=0 



(2.20) 



is conserved. Notice again that the charge makes sense provided the defect is not too 
severe with \u — v\ x=0 < a. 
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Notice that adding and subtracting the two defect conditions (2.11) gives the 
following pair of relations at x = 0: 

u t -v t = i (u x + v x ) Q + i(u - v) (\u\ 2 + \vf} 

u x — v x —(u + v)Q. (2-21) 

If a pair of conditions such as these held for every x, and not merely at x — 0, (2.21) 
would be recognised as a Backlund transformation: differentiating the second equation 
with respect to x, using the second equation to substitute for u x —v x wherever it appears, 
and adding i times the first equation gives 

iu t + u xx + 2w|-u| 2 = iv t + v xx + 2v\v\ 2 . (2.22) 

On the other hand, differentiating the first with respect to x and the second with respect 
to t and subtracting leads to 

iu t + u xx + 2-u|w| 2 = — (it>t + v xx + 2v \v | 2 ). (2.23) 

Hence, combining these manipulations, both u and v satisfy the NLS equation. In fact, 
this is exactly the Backlund transformation for NLS given by Lamb [17] (see also [18], 
[19]). 

3. Argument for integrability 

In order to verify that introducing a jump-defect does not destroy integrability it 
is necessary either to investigate additional conserved charges or, more effciently, to 
examine a suitably modified Lax pair. The latter approach is generally superior since the 
Lax pair provides a generating function (as a Laurent series in the spectral parameter) 
for an infinite set of independent conserved quantities. 

A Lax pair for the NLS has been provided in [23] (see also [20]). For the above 
choice of conventions, a satisfactory Lax pair is: 

L{u) = i (ucr+ + ua ) + A03 

M(u) = i(2uu + A 2 ) 03 + (u x — \u)a+ — (u x + \u) <j_ , (3.1) 



where 



with the property 

d t L - d x M + [M, L] = \u t + u xx + 2w|w| 2 = . (3.3) 

Moreover, (3.3) holds independently of the choice of spectral parameter A. 

Following the ideas introduced in [8] and [22], modified Lax pairs can be devised 
which will build in automatically both the bulk equations and the defect conditions. In 
designing these it is natural to introduce two extra points a < and b > which are 
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the endpoints of two regions overlapping the defect, one on the left R , — oo < x < b, 
and one on the right R + , a < x < oo Then, a suitable pair can be defined as follows: 

L~ =L(u)6(a-x) 

( T ; O (n. _ ni\ 

2 + M 2 ) 



M 



M(u) + 9(x - a) 



i (-u — -u) 



2ft 



2Q 



u 



(7_ 



2ft 



(u t - v t ) - —{u + v) + 



It 



2fi 



■d«l 



<T 4 



(3.4) 



and 

L + 

M + 



L{v) 6{x - b) 
M{v)+6{b- 



v 



i I \ ^ i \ ( M ~ V ) n i2 



2fi 

+ M 2 ) 



(i«r+i^i 2 ) 



(7_ 



(3.5) 



These modified Lax pairs provide the equations of motion together with all the defect 
relations as a consequence of zero curvature conditions of the type (3.3). In the 
overlapping interval a < x < b, the two matrices M + and M~ must be x-independent in 
order to maintain the zero curvature condition, though not necessarily equal. Rather, 
they must be related by the following 'gauge' transformation (see [22] for more details), 

d t )C = )CM + (b,t)-M-(a,t))C. (3.6) 

Note that since M ± are ^-independent this implies that the field u and v are also x- 
independent in the overlapping interval. It can be verified that the following choice for 
K 



1 



/C = /+- 
A 



-\{u — v) 

-n 



(3.7) 



i(u — v) 

where I is the 2x2 identity matrix, works perfectly, in the sense that (3.6) is identically 
satisfied. In the limit a — > the discontinuity tends to zero and /C — > /, as one would 
expect. 

If one was to consider instead the conditions associated with a 5-impurity, namely, 

u x = v x — (tv, v x — u x + era, x — 0, (3.8) 

a similar construction would lead to the conclusion that there was no suitable matrix K 
and hence no Lax pair for this type of defect. This suggests the 5-impurity is not 
integrable, a fact apparently consistent with numerical studies [5], [6]. One could 
envisage more general linear defect conditions, but the conclusion remains the same. 

Once /C is given, and assuming suitably decaying fields at ±oo, or periodic boundary 
conditions, conserved quantities will be generated by 



Q(A) = Tr 



Pexp(^J dx L(u)^j /CPexp^y dx L(v 



(3.9) 



The final and necessary step would be to demonstrate that the charges generated by 
(3.9) are independent and in involution. However, the discussion of this topic will be 
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- u x u) t - ^{uu t - UtU) 3 



(2uuA + iu x B — iu x C) x — (2\u u + iu xx )B + (2\u u + \u xx )C . (3.13) 



postponed. In the meantime, following Ablowitz et al. [24], it will be shown how the 
Lax pair can be used to generate conserved charges using the integrability constraints 
represented by the zero curvature condition. 

The bulk Lax pair for NLS can be written alternatively as follows 

where the entries of the matrix M can be obtained from (3.1). Applying the zero 
curvature condition (3.3) to this Lax pair, the following relations are obtained 

A x = iuB -iuC, XB = 2mA - B x + iu t , AC = 2\uA + C X - \u t . (3.11) 

These allow the conserved charges to be calculated as the coefficients in an expansion in 
1/A. To see how the argument proceeds, take a couple of steps in detail. Substituting 
the last two equations of (3.11) into the first of (3.11) it follows that 

A x = -- [(uu) t + i(uB + uC) x - i{u x B + u x C)} . (3.12) 
A 

Repeating the substitution from (3.11) one finds 

A x = - j [{uu) t + i(uB + uC) x \ - ^ 
1 

The process can be repeated to obtain further elements of the expansion in higher order 
in 1/A. However, (3.13) is enough to illustrate the point. If total ^-derivatives integrate 
to zero (assuming u and its derivatives vanish as \x\ — > oo) then integrating (3.13) over 
— oo < x < oo and concentrating on the first two terms gives 
d f°° _ d f°° 1 

— / dx(uu)=0, — / da; -{uu x — u x u) = 0. (3.14) 
dtJ-oo v ; dtJ-oo 2 V ; K J 

These are the 'probability or number' and momentum densities, respectively, used 

previously. 

With a jump-defect the situation is different and it becomes necessary to deal with 
the modified Lax pairs defined in two regions whose overlap contains the defect. In the 
matrix form, the Lax pair to use reads as follows 

where the entries of the matrices M + and M + will not be specified here but they can be 
obtained from the formulae (3.4) and (3.5). The Lax pairs (3.15) and (3.16) are defined 
in the regions R~ and R + , respectively. Within the overlap the zero curvature condition 
requires 

d x M~ = d x M + = 0, (3.17) 
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and therefore is constant throughout the overlap although the two constant values 
are not necessarily the same. In fact, the 'gauge' transformation (3.6) supplies the link 
between the values of these constants via 



M + (b,t) =JC~ 1 d t lC + JC- 1 M~(a,t)IC, 



(3.18) 



Next, imagine shrinking the overlap region by allowing its endpoints a, b to approach 
the defect location at x — 0. In these circumstances, the strategy used before, to obtain 
the conserved charges in the bulk, can be applied. For example, the analogue of (3.12) 
is 



A x + A+ = -- [{uu) t + (vv) t + i(uB x + uC x ) x + i(uB+ + uC+) 3 

-i(u x B~ + u x C~) - \{u x B+ + u x C+) 
Integrating over — oo < x < oo gives 

1 



(3.19) 



1 


-g ,0 


x=o ~ A 


dt J-oo 



d_ 
di 



dx (vv) 



A 



i(uB~ + uC~) - iuiK^dtK. + Kr x M~lC) 12 - m{Kr x d t K + Kr x M-K)n 



1 



/0 „ poo 

dx\{u x B~ + u x C~) + / dx\(u x B+ + u x C x ) 
-oo JO 



x=0 

(3.20) 



where the elements of the matrix (3.18) appear. In turn, these can be obtained using 
(3.4) and (3.7). Doing so, the first two lines of the integrated expansion (3.20) can be 
calculated and, schematically, the result is the following 



i r .9 

x=0 + A I dt 



dx (uu) + / dx (vv)) + 
Jo J 



J x=0. 



1 

A 



/0 * poo 

dxi(u x B~ + u x C x ) + / dxi(u x B+ + u x C^ 
-oo JO 



1 

0, 



\x=0 



(3.21) 



where the unwritten parts contain terms in u, v and their spatial derivatives. These 
terms cancel out after a lengthy calculation. In the first line of (3.21) it is possible 
to recognize the correction due to the defect, which has to be introduced to make the 
'number' charge conserved (2.19). The iteration process has been checked up to order 
1/A 2 confirming the correction to the momentum already indicated in (2.12). The details 
of this verification are straightforward and will be omitted. 

It should be clear that these arguments can be repeated for an arbitrary number of 
jump-defects placed along the a;-axis. Each defect is treated locally, introduces its own 
parameter, and further complicates the generating functional (3.9) for the conserved 
quantities. 



4. A single soliton meeting a defect 



It is interesting to investigate what happens as a soliton approaches a jump-defect. 
With the above normalisation for the nonlinear term in the field equation, a one-soliton 
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solution is given by 
2o EF 

u= — , E = exp [a(x — 2ct)\ , F = exp [i (cx + (a 2 - c 2 ) t)] (4.1) 

1 + E 2 L V V / / J 

where a > and c are free, real parameters. For general values of these parameters, the 
modulus and phase of the soliton travel with different speeds. A slightly more general 
solution is obtained by shifting the origin of x, for example, and multiplying u by an 
arbitrary phase. 

When there is a defect, the field v on the other side of it is taken to be given by a 
similar expression but with different shifts in the modulus and the phase components; 
in other words, take v to be 
2apqEF 

V ~l+p^E^ [ ^ Z) 
In view of the form of the defect conditions it is also reasonable to suppose p is real. 
Then, the first miracle which should occur concerns the square root in the defect 
conditions whose argument will have to be a perfect square. For this one already requires 
the following two relations between p and q: 

|1 — pq\ A a 2 
[1-p 2 ) 2 = ~a~ 2 

Then, checking the second of equations (2.21) requires 



a + ic — a 

pq = —. - , 4.4 
a + ic + a 

and therefore 

9 (a — a) 2 + c 2 a — a + ic a + a — ic , 4 . 

V = , , 2 ^ 2 , 1 = | r | ~ , ■ I • (4-5) 

(a + Qt) z + c z \a — a + ic| |a + a + ic| 

The latter is quite a nice expression for p since it can never vanish (for c ^ 0), and it 
approaches unity as c — > oo. Clearly, the soliton cannot be 'eaten' by the defect and 
the faster it travels the less it is affected. Also, p is not sensitive to the sign of c, but 
changing the sign of c replaces q by q. As \a\ — > oo the parameter g — > — 1, indicating 
that when the defect parameter is sufficiently large the soliton will have inverted its 
shape when it emerges from the defect (but, since p — > 1, the inverted soliton will not 
be significantly delayed). As a — > 0, both p, q — > 1 and the effect of the discontinuity 
disappears, as it should. Notice, the above collection of formulae refer to the positive 
square root Q; changing the sign of the square root will require changing a — > —a in 
the above expressions. 

Overall, the picture for NLS is quite similar to that discovered previously for sine- 
Gordon except that in the sine-Gordon case a soliton may be absorbed by the defect, 
or converted to an anti-soliton [7] according to the choice of defect parameter. These 
possibilities cannot occur for NLS. 
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5. The two-soliton solution and a defect 

One would expect that a pair of initially widely separated solitons approaching the 
defect should pass through it independently, each experiencing the jump-defect as if the 
other soliton were not there. Indeed this is the case although to verify it explicitly is a 
formidable calculation requiring the use of Mathematica or Maple. In this section the 
principal steps in this verification will be outlined. 

The multi-soliton solution has been given in closed form by Zakharov and Shabat 
[23] (see also [20]) in terms of the scattering data of the linear problem associated with 
the NLS. This data is a set of complex constants Aj = Cj + idj dj > and functions 
Tj(x,t) = 7j expi(AjX — \ 2 t) where jj provides the initial position and phase of the 
soliton labelled j. The two-soliton solution is then described as follows 



u 
2 

where 









aiAiri + a 2 A 2 r 2 


: + aiAiTi 


r 2 i 


2 + a 2 A 2 T 2 


ri| 


2 


1 + 




2 


r 2 ; 


2 + 


Ail 


2 




2 + 


a 2 ; 


2 


r 2 i 


2 _ 


4a 1 a 2 (r 1 r 2 + r 1 r 2 )/|A 1 -A 2 | 2 



(5.1) 



_ (A 2 - Ai) _ (Ai - A 2 ) , v 

A2 " (A^AO ' Al " (aT^) • (5 ' 2) 
Setting c — {c\ — c 2 ), a solution of this type on the left of the defect can be written 
conveniently as: 

u = a 1 ^ 1 F 1 [A + (l + E 2 2 ) - 2ica 2 (l - E 2 2 )} - a 2 £ 2 F 2 [A_(l + g) - 21^(1 - gfj] 

2 ^ + 5l)(l + ElEl) + ^ + 5l)(El + El)-Aa 1 a 2 E l E 2 (F^ + Fl)/F l F 2 [ ' } 

where 5 + = (ai + a 2 ), 5_ = (ai — a 2 ), A± = 5 + 5- ± c 2 and 

Ej = exp [aj(x — Cjt) — aj x j], Fj = exp [i(c,- x + (a 2 — c 2 )t + O j)] j = 1, 2. (5.4) 

Note that expression (5.3) contains the one-soliton solution (4.1) on setting, for instance, 
a 2 = c 2 = 0. 

To the right of the defect the expression (5.3) will be modified as follows 
v a 1 C 1 [A + (l + H 2 2 ) - 2ica 2 (l - Hp] - a 2 G 2 [A_(l + Hi) - 2\ca l {l - Hf)] 

2 (c 2 + 5i )(i + mm) + (c 2 + 5 2 )(^ 2 + m) -^Moim + cim)/^^) 1 ■ ' 

where Gj = ZjEjFj, H? = WjEj and Wj = p 2 , Zj = pj qj with j — 1,2 and \q±\ = \q 2 \ = 1. 
The constants p^, represent the delays in position and phase, respectively, for the two 
solitons residing inside the solution v, by analogy with the notation used for the one- 
soliton solution on each side of the defect. 

The first step in exploring the consequences of the defect conditions is to find the 
circumstances under which the argument of the square root f2 = (a 2 — \u — v\ 2 ) 1 / 2 is a 
perfect square. The most general polynomial whose square can match the argument of 
the square root must have the following form, 

n = {a + 61 El + b 2 El + di E{ + d 2 E\ + ei E\El + e 2 E 2 X E\ 

+ 9l E\E 2 + g 2 E X E\ + hE 1 E 2 + r E\El + s E{E\ + / E 2 X E 2 2 ). (5.6) 
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The coefficients of Q 2 are vastly overdetermined and it is expected that just a few 
relationships among the parameters will suffice. Indeed, this turns out to be the case, 
and the following constraints on the constants Wj, Zj are all that is required, 

[ a .( z . _ !)(„,. _ Zj )f = [a(l- Wj ) z 3 \ 2 j = 1,2 (5.7) 

(ai + a 2 -\c)(wi- zi)(z 2 -l) =2a(zi-wiz 2 ) a* ^ 0, i — 1,2 
(ai — a 2 + ic)(zi — l)(z2 — 1) = 2a(zi — z 2 ) a, ^ 0, i — 1,2. (5.8) 

Note that, contrary to relations (5.7), the expressions (5.8) are not real and therefore 
their complex conjugate partners have also to be taken into account 

(ai + a 2 + ic)(z 1 - l)(w 2 - z 2 ) — 2a(z 2 - w 2 z^) a { ^ 0, i — 1,2 

(ai - a 2 - ic)(iui - z 1 )(w 2 - z 2 ) — 2a(w 1 z 2 - w 2 Zi) a, 7^ 0, i — 1,2. (5.9) 

Implementing these constraints, the second part of the defect condition (2.21) is satisfied 
provided the following two relations hold 

2* = — J = 1,2 5.10 

a j + a + i Cj 

and therefore 

(a,- - a) 2 + c 2 

w 3 = )\ ' 3 2 J = 1,2. 5.11 
(a,- + a) 2 + cj 

Matching the conditions in this way demonstrates that solitons are transmitted through 
the defect independently of one another. It is expected this property should hold for 
arbitrary numbers of solitons although there is no general proof of that yet. 

In the sine-Gordon model it has been remarked that the delay experienced by a 
soliton passing a defect is actually the square root of the delay that would be experienced 
by the same soliton passing another whose rapidity was equal to the defect parameter 
[13]. Using the following change of variables suggested in [23] (see also [20]) when 
ci > c 2 , 

t+ = t 1 x 1 17 = ^ r+ = ^ r 2 - = r 2 A 2 , (5.12) 

Ai A 2 

the two soliton solution (5.1) may be rewritten as follows 



u 
2 









+ a 2 r 2 




f a 1 T 1 


r 2 | 


2 + o 2 r+i 


rfl 
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i + |r+| 


2 + 




2 + 




2 


1 2 


2 _ 


8a 1 a 2 Re[(r+r 2 )/(A 1 -A 2 ) 2 ] 



(5.13) 



Note that the change of variables (5.12) always holds provided Ai ^ A 2 , A 2 . Moreover 

ir+pir^i 2 = |rr| 2 |r 2 1 2 = |r x | 2 |r 2 1 2 - (5.14) 

Functions Tf(x,t) are defined as follows 

Ff(x, t) = ~ff exp [\(\jX - \ 2 t)}, ~ff = exp (a jX ± + i0± ) (5.15) 

where now 

Ap = exp [-ai(x^ - x m )} = exp [a 2 (x£ 2 - ^02)], 

Ag i = exp [-i(0^ - O1 )], Ag 02 = exp [i(<f>% 2 - (j) 02 )} (5.16) 
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A ^ = ^—T^——- 2 . (5.19) 



are the two-soliton scattering data. In fact, looking at the limit t = ±00 of the two- 
soliton solution (5.13), for instance along the first soliton trajectory x — C\t — constant, 
it is possible to obtain the following one-soliton solutions (ci > c 2 ) 

u\ t-,00 - 1 - |p+p u\ ^_oo - - - 2 . [b.l <) 

Using (5.12), explicit expressions for the scattering data Ap Q , Ag i and Aq 02 are: 
a (a 1 -a 2 ) 2 + c 2 

A P° = 7 — 7 — yT~T^2^ 5 - 18 
(ai + a 2 y + c 2 

\-i ,+ - a? — a, + 2ia 2 c + c 2 

(Ag 01 )- = exp [i(0+ - 00,)] = - _ — - - , 

a 2 — a\ + 2iaic + c 
a i — a i ~ 2iaic + c 

Note that (Ag i) _1 can be also written as follows 

! _ (a 2 -a2)2 + c 4 + 2c 2 a 2_ 6c 2 a 2 +4ia2 ( c 2 + a 2_ a 2) 

(Agoi) - {al-alY + c± + 2c\al + al) " (5 ' 20) 

It is then easy to check that Ap and (Ag i) _1 coincide with the squares of the shifts 
in the modulus p and the phase q (4.5) experienced by the soliton progressing through 
a defect, provided c 2 = and a 2 = a, with a being the defect parameter. 

6. Bound states 

With a delta function impurity, in a free field situation with a suitable coupling, there 
will be a bound state. Curiously, this is not the case for the free field limit of the sine- 
Gordon model with the jump-defect conditions provided by (1.2), nor is it the case for 
the full sine-Gordon model with a jump-defect. In the sine-Gordon model this type of 
purely transmitting defect does not permit bound states for any value of the coupling a. 
Allowing a jump-defect in the quantum sine-Gordon model does have interesting effects, 
however, and it appears a defect can be excited, albeit with a finite decay width. This 
is the quantum analogue of the classical property permitting a jump-defect to swallow 
a soliton [13]. 

On the other hand, a jump-defect in the NLS model certainly has bound states 
associated with it. The first indication of this arises using the linear, potential-free, 
Schrodinger equations and the linearised versions of the defect conditions (2.21) 

i (x 
u x = - —{u t - v t ) + -(u + v) 
la I 

i ol 

v x = - —{u t - v t ) - -{u + v). (6.1) 
It is easy to check these allow the travelling wave solution 

u = u exp (— \k 2 t + ikx), v = v exp (— \k 2 t + ikx), v — 7 — u , (6.2) 

k — ia 
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where k is real, from which the existence of bound states associated with the defect may 
be deduced. If a > (< 0), (6.2) indicates a bound state when k = ia (— ia) for which 
there are the square integrable solutions 

u — 0, v — vq exp (ia 2 t — ax) ; (u — uq exp (ia 2 t + ax), v — 0), (6.3) 
respectively. 

Interestingly, these bound states have their counterparts in the fully nonlinear 
version of the model. For example, taking a to be positive, with a — a, c = in 
(4.1), the field configurations 

2ae a(x-x ) Q ia 2 t 

u = — —. -, v = 0, 6.4 

I _|_ e 2a(x~x ) ' v ' 

in the regions x < and x > 0, respectively, satisfy the jump defect conditions. If the 
parameter xq = 0, the quantity Q actually vanishes at the defect, as does the momentum 
of course, while N = a and E = — a 3 /3 are exactly half the values the number and 
energy would have had for a static breather solution (with a — a, c = 0) in the bulk. 
In fact, the latter values for iV and E remain correct if x ^ 0; in those cases f2 does 
not vanish but compensates the bulk integrals. Note, this solution is entirely consistent 
with the expressions discovered before in (4.5), noting that if c = then p vanishes 
precisely when a = a. 

7. The modified Korteweg-de Vries 

One might wonder if the modified KdV (mKdV) equation, or the KdV itself, both of 
which certainly possesses Backlund transformations, also allow discontinuous solutions, 
or defects, provided suitable conditions are imposed. 
In the bulk, a mKdV equation is [21] 

v t - 6v 2 v x - v xxx = 0, (7.1) 

or, setting v — p x , 

Pxt ~ GPIPXX ~ PXXXX = 0, (7.2) 

the latter being suitable for a Lagrangian description with Lagrangian density 

£ = IjPxPt - ^(Px) 4 + ^(Pxx) 2 - (7.3) 

Integrating (7.2) once with respect to x, and assuming all derivatives are asymptotically 
vanishing gives an alternative: 

p t - 2pl - p xxx = 0. (7.4) 

There is another version of mKdV (obtained by the replacement v — > if) in which 
the middle term on the left hand side of (7.1) changes sign. However, it is equation 
(7.1) which allows real soliton solutions (see below), not the alternative version [24]. 

Over the whole line, the quantity 

I roc 

P = - dxp 2 x , (7.5) 

Z J—OD 



Jump-defects in non-relativistic field theories 



16 



is conserved (with the usual assumptions at ±00) because 

where the last expression made use of the alternative equation (7.4). 

Next, suppose there is a 'defect' at x — 0, with fields p, q on either side of it. The 
quantity P defined by 

1 r0 1 roo 

P= 2/00^+2! dXq * ^ 

will not be conserved but, as a consequence of (7.6), its time derivative will be related 
to a boundary term as follows, 

«=(^-T-fL-(«--!-fL- ™ 

The question is how to write the latter as a time derivative of the fields and their 
derivatives evaluated at x — 0. 

To see where the defect conditions are coming from, it is necessary to return and 
consider the complete action 

^ = ydt|y° dx£(p) +B + J™dxC(q)}, (7.9) 
and its variation with respect to p or q. Thus, for example, varying p gives, 

6A= J dt^J dx (^bptPx + ^SpxPt - 25p x pl + $PxxPxx^J (7.10) 

/. dB dB . 0B\ ) 
+ (5p t — + 5p— + 5p x —) , 7.11 
V dp t dp dpxj x=0 ) 

which, on integrating the first term by parts with respect to t and x, keeping the 
boundary terms in x, and setting the variation to zero, leads to the field equation for p 
together with 







' (I n 3 dB d dB\ ( dB 

OP «Pt ~ 2 Px ~ Pxxx H 7\ 7777^— \+0Px\ Pxx + 



(7.12) 

x=0 



op at opt ) \ op^ 

At the defect, there is no necessity for bp and 8p x to be related but (7.4) can be used 
(in a limiting sense) to eliminate p XX x- Hence, the defect conditions on p should read 
1 dB d dB dB . . 

Similarly, the relations to be satisfied by q at the defect are: 

n 1 dB d dB n dB fr7iA . 

° = 2* + ^-aaS' ° = - q - + dq-x- (7 - 14) 

The next step is to find a suitable defect term B. However, this is not quite 
straightforward and before doing so it is worth making a remark. 

It is tempting to assume the defect term depends only on p, q and their first 
derivatives (in fact this was tacitly assumed above). However, there is no reason in 
principle why this term should not depend on p xx and q xx since the field equations 
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cannot at the defect relate these derivatives to anything else (but note, the same would 
not be true of p xxx , or q xxx ). Moreover, since the Lagrangian (7.3) depends upon both 
first and second derivatives, it should not be surprising that defect conditions should 
also involve higher spatial derivatives. If B does depend on these derivatives then 
there will be two more defect conditions coming from varying the boundary term alone, 
unbalanced by anything in the bulk; they are 

|^ = 0, |*=0. (7,5) 

OPxx OQxx 

On the other hand, a Backlund transformation for the mKdV equation was found 
long ago [17, 25] and can be written in the following symmetrical manner: 

Px + q x = asin(p - q), 

Pt + q t =ol (p xx - q xx ) cos(p - q) + a (p 2 x + qfj sin(p - q) 

= ct (p xx - q xx ) cos(p -q) - 2ap x q x sin(p - q) + a 3 sin 3 (p - q). (7.16) 

The parameter a is arbitrary and the pair of equations is symmetrical under 
interchanging p and q, and simultaneously making the replacement a — > —a. In the 
bulk, the derivative of the first of the pair (7.16) is 

Pxx + qxx = ct (p x - q x ) cos(p - q), (7.17) 

which will provide the extra equation coming from (7.15) (see below). Note that, under 
the circumstances being explored here, the ^-derivatives are frozen, therefore (7.17) 
cannot be a consequence of the first of conditions (7.16). 

Using (7.16) and (7.17) the expression (7.8) can be simplified to 

Pt = j t [-<*cos(p-q)] x=0 , (7.18) 

and therefore, 

P + «[cos(p - q) - 1} X=0 (7.19) 

is conserved. The constant has been chosen to ensure the momentum stored at the 
defect is zero when there is no discontinuity. It is worth noting that because of the 
discontinuity the charge that in the bulk is obtained from the density Af = p x is not 
conserved. The same could be said for the KdV model discussed in section 9. 
As previously, it is useful to put 

B = ^(qp t -pq t )-B (7.20) 

and the defect conditions may be rewritten 

dB dB OB dB dB dB ,„ , 

(Pt + Qt) = -2^- = 2— ; p xx = — , q xx = -— ; = - — = - — . (7.21) 

dp dq dp x dq x dp xx dq xx 

Putting everything together, a suitable choice of B appears to be to take 

B = - {p xx - q xx ) (p x + q x -ct sin(p - q)) 

+ ^ cos(p - q) \p 2 x + q 2 x - 4p x q x + a (p x + q x ) sin(p - q) + a 2 ] . (7.22) 
6 L J 
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Using (7.22) the defect conditions are exactly equivalent to what would be the Backlund 
transformation in the bulk together with (7.17) were the conditions not frozen at x — 0. 
Finally, it is worth remarking that the mKdV defect potential is not a simple function 
of the two fields u and v and the parameter a but rather of the 'potentials' p and q. 



8. The mKdV soliton passing through a defect 

The single soliton for the mKdV equation can be conveniently written in terms of p in 
the form 

1 + \E 

e ip = — , E = exp \a(x - x + a 2 t)} , (8.1) 

1 — iE <- J 

where a and exp (x ) are both real parameters. With the conventions adopted in the 

last section, the soliton moves from right to left along the x-axis. Since the derivative 

of the above expression is either always positive (E > 0) or always negative (E < 0), 

there is also a matching anti-soliton obtained by replacing E — > — E in (8.1). This is 

easily seen, on noting 

2aE . . 

» = TTEi- (8 ' 2) 

If there is a defect, the soliton on the other side of it will have the form 

eig = l + lzE ? E = exp \a(x - x + a 2 t)] , (8.3) 
1 — izE l J 

where z is a parameter to be determined by the defect condition. 

The first of the defect conditions readily reveals that 

z= . (8.4) 

a + a 

This appears to suggest that a slow soliton is not affected much whereas a fast soliton 
has z < 0. This means that a fast soliton flips to an anti-soliton (E — > —E in the 
expression (8.1)). When a — a, the soliton is eaten. This behaviour is very similar to 
that of a soliton meeting a defect in the sine-Gordon model where all these effects are 
similarly apparent. Despite the non-locality, there appears to be nothing particularly 
pathological about this case. As a final remark in this section, it is clear the effect of 
the defect disappears as a — > oo. 



9. The Korteweg-de Vries equation 

It is not difficult to repeat these steps for the KdV equation. However, there does 
seem to be some curiously different behaviour that will become apparent as this section 
proceeds. 

In the bulk, the KdV equation is [21] 

u t - 6uu x + u xxx = 0, (9.1) 
or, setting u = p x , 

Pxt ~ QPxPxx + Pxxxx = 0, (9.2) 
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the latter being suitable for a Lagrangian description with Lagrangian density 

C=\pxPt-(Px?-\(Pxx) 2 - (9.3) 

Integrating (9.2) once with respect to x, and assuming all derivatives are asymptotically 
vanishing gives an alternative: 

Pt-3p 2 x +p xxx = 0. (9.4) 

Over the whole line, the quantity 

I poo 

P = - dx(p x ) 2 , (9.5) 



2 

is conserved (with the usual assumptions at ±oo) because 

(§) t = { 2pl - pxpxxx + *f) = { ptpx - pi+v -f) ' (9 - 6) 

where the last expression made use of the alternative equation (9.4). 

Next, suppose there is a 'defect' at x = 0, with fields p, q on either side of it. The 
quantity P defined by 

P=~\ dxpl + - dxq 2 x (9.7) 

Z J-oo Z JO 

is not conserved but, as a consequence of (9.6), its time derivative will be related to a 
boundary term as follows, 

Pt=(ptp x -Pl + P f) -(ft&-<£ + ^) ■ (9-8) 

The question is how to write the latter as a time derivative of the fields or their 
derivatives evaluated at x — 0. 

The right hand side of (9.8) can be organised a little differently (dropping the 
explicit reference to the point x — 0, which is to be understood from now on) to, 

-(Px - qx){pl + p x q x + ql) + ^{Pxx ~ q X x){Pxx + Qxx) + PtPx ~ QtQx, (9.9) 
and then simplified by setting 

Pxx + Qxx = (p - q)(Px - q x ) 

P 2 X + PxQx + q 2 x = ~[Pt + qt + (P~ q)(Pxx - q X x)] , (9-10) 



to get 

1 



2 (Pt-qt)(Px + qx)- (9.11) 



Finally, setting in addition 



p x + qx = 2a + ^(p-q) 2 , (9.12) 



one finds 

P- d 



«(p - q) + j^ip - q) 3 



x=0 



(9.13) 
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where the parameter a is arbitrary. The first of the conditions (9.10) would follow from 
(9.12) in the bulk but, here, as already mentioned before, the x-derivatives are frozen. 
The second of equations (9.10) together with (9.12) provides a Backlund transformation 
for KdV in the bulk; at least in the sense that if p satisfies (9.2) so does q (or vice- 
versa), for any choice of a. This is the form of Backlund transformation constructed 
by Wahlquist and Estabrook [26] for generating multi-soliton solutions to the KdV 
equation. 

This time using the Lagrangian together with a defect contribution (assuming the 
latter depends on p and q and their derivatives p x , q x and p xx , q xx ), leads naturally to 
the defect conditions: 

1 <9B d dB dB dB 

° = °=- p " + w,' 0= as:- (9 ' 14) 

Similarly, the relations to be satisfied by q at the defect are: 

ft 1 dB d dB n dB n dB 

0= o^ + TT-^T^-' = fc + — , = -—. 9.15 
2 dq dtdq t dq x dq xx 

A suitable choice for the jump-defect potential might be the following 

B = j (qPt - pqt) + (pi + p x q x + ql) {p-q) + \ (Pxx ~ q xx ) 



p x + q x -2a- ^{p-q) 2 



- 3 (p x + q x ) (p -q)(a + ^(p - qf) + Ga 2 (p - q) + 2a(p - q) 3 + — (p - qf. (9.16) 

This looks quite complicated but it does not seem easy to find anything simpler that 

would be able to provide the relations needed. In this respect, both the last equations 

in (9.14), (9.15) give precisely the equation (9.12), while two of the other conditions, 

still from (9.14), (9.15), give 

1 3 
Pxx= ^(Pxx-qxx) + (2p x + q x )(p~q)-3a(p-q)--(p-q) 3 , (9.17) 

1 3 
qxx = - - (Pxx - qxx) - (Px + 2q x ) (p-q) + Sa(p - q) + -(p - q) 3 . (9.18) 

Adding these gives 

Pxx + qxx = (px-q x )(p-q), (9.19) 

which, in the bulk, is just the derivative of (9.12); subtracting them gives (9.12) again 
(the second derivatives exactly cancelling out). 

Finally, the other pair of equations from (9.14), (9.15) give effectively the same 
condition, namely 

2 (Pt + (It) = (Pxx ~ q xx ) {p-q) + (p 2 x + p x q x + ql) ■ (9.20) 

This is not quite straightforward to see and makes use of (9.12) again. Indeed, the 
coefficients of the polynomial part of B were chosen to ensure this worked out. 

With this particular choice of B, the 'momentum' is compensated in precisely the 
manner envisaged before (in (9.13)), and 

1 



a(p - q) + -^{p - q) 3 



(9.21) 

x=0 
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is conserved. Moreover, the parameter a is completely free. 

The next bulk conserved quantity, the 'energy', has a density 

{Pxf + \{Pxxf 

It is not conserved but the defect contributes exactly 
provided the defect potential (9.16) is written 

B = \{m -pqt) - B. 

Hence, using the boundary conditions 



E + B = E 



(pl + Pxq x + q 2 x ) (p-q) + (p- q) 3 (a + — (p - g) 2 ) 



x=0 



(9.22) 



(9.23) 



(9.24) 



(9.25) 



is conserved. 

The situation is fairly similar to what was found for mKdV. It is curious that 
neither of the compensating terms can be rewritten in terms of the the original KdV 
variables u on the left, or v on the right (u — p x , v — q x ). The expressions for energy and 
momentum reveal that there is nothing particularly special about a vanishing parameter 
a. 

In general, provided the defect potential B is chosen as it is in (9.24), the momentum 
is conserved provided 



'dB_ 

dp x 



- ( ^ ) +(p*-qx) 



for mKdV and 



' dB 
dPx 



dB 
dq. 



-[-*—) -(Px-Qx) 



dB 

2^ + (pl + q 2 x)(Px + qx) 



2— + 2(pl+p x q x + q 2 x ) 



0, (9.26) 



0, (9.27) 



for KdV, respectively, where 

2 



OB 

dp x 



dB 

dq x 



ipxx ~\~ qxx ) ipxx qxx ) j 2 



dB 

dp 



together with 



(Pt ~ qt){Px + q x ) = 2^, 
for either model, where T is a function to be determined. 



ip t + qt), (9.28) 



(9.29) 
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10. KdV solitons and the defect 

A single soliton can be described conveniently by choosing 
2 clE 

p = Po — — , E = exp a(x — x — a 2 t) = p exp a(x — a 2 t) , (10-1) 

1 + E L J L J 

where po^^o are constants. In terms of the original variable u the solution has the 
well-known characteristic 'bell-shaped' form 

2a 2 E , 

« = * = -(TW < ia2 > 

The parameter a is real and the soliton is independent of the sign of a (since p x is the 
same whatever the sign). Inevitably, the soliton moves to the right. 
Then, it is not difficult to check that picking a similar form for q 

2 a <j E r . \i /\ 

q = g - ^ — S> ^ = ex P ^ ~~ a *) ' ( 10 - 3 ) 

1 + CT-C/ L J 

and redefining the defect parameter a = — /9 2 /4, equation (9.12) implies 

At least that is the case assuming the positive square root is taken for po — qo- The 
other equation (the second of (9.10)) is just an identity using only the fact 

Po - Qo = a \^— \ , (10.5) 
\p-o-J 

as can be checked easily using an algebraic computing package. Finally, the first of 
conditions (9.10) is also an identity and can be checked directly. The relation (10.5) 
does not depend on taking square roots. 

This implies a soliton encountering the defect, provided it is not travelling too 
quickly, will be delayed. When a approaches f3, a tends to zero and the soliton will 
have been 'eaten' (since q = qo). There is a mystery if the soliton is moving too quickly 
because the solution on the right of the defect develops a singularity. This is puzzling 
for another reason. If the defect parameter is taken to be very small then only a very 
slow soliton will remain non-singular, meaning that in the limit as the parameter goes 
to zero one cannot recover a single bulk soliton because the speed of the non-singular 
soliton would also have to approach zero. Similarly, in the limit of large defect parameter 
(Po — Qo) 2 00, implying that the effect of the defect on the soliton does not entirely 
disappear, although the fields to either side of it will be the same (because p x — > q x ). 
In these respects, the KdV defect behaves curiously differently from those encountered 
elsewhere. 

One might wonder whether having a single soliton on the left of the defect and 
a double soliton on the right might provide some hints towards solving the mystery. 
To explore this, consider the following single one-solution on the left and two-soliton 
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solution on the right of the defect, respectively 

2 CL2P2E2 



P = Po 



1 + p 2 E 2 ' 



= + (oj ~ ^2) (1 + Pi^i + ff P2^2 + o-pip 2 E 1 E 2 ) , 1Q g . 

° (ai - a 2 )(l - apip 2 E 1 E 2 ) - (ai + a 2 ){p\E 1 - ap 2 E 2 ) ' 



with 



PjSj = exp aj(x — x j — a 2 t) j — 1,2. (10-7) 

As pointed out in [26], a two-soliton solution is non-singular provided one of the two 
component solitons is actually singular. For instance, in (10.6) suppose a 2 > cti, then 
for q to be regular pi > and op2 < and consequently the faster soliton is singular. 
Checking the Backhand equation (9.12) for the solution (10.6) reveals the following 

al = (3 2 , (p -qo) 2 = (3 2 , a = 1; (10.8) 

therefore, provided the incoming soliton is regular (p 2 > 0), the resulting two-soliton 
solution is singular. One could think that the 'fast' soliton is trapped by the defect since 
it does not have enough energy to escape. On the other hand, if the one-soliton solution 
on the left is singular (p 2 > 0) then the resulting two-soliton solution would be regular 
(still with a 2 > ai). The singularity of the incoming soliton can be kept to the right of 
the defect by a suitable choice of the constant X02, in fact, simply taking x 02 > will 
suffice. Such a singularity remains on the right of the defect as t — > 00. 

The KdV equation also allows other progressing singular solutions of the following 
kind, 

2 c 2 c / _ 1 \ , . 

u = — — -, p — [x — xi — -ct , 10.9 

(x - x l - ct) 2 6 F x-Xi-ct 6 V 27' v ; 

where c, x±,Xi are constants. The jump-defect will also affect these in the following 
manner. Let q be given by 

-(x-x 2 --ct) , (10.10) 



x — x 2 — ct 6 V 2 
where x 2 ,x 2 are two additional constants. Then the defect conditions require 

xx- 12 ( ° + 1 \ (10 11) 

c{xi—x 2 y \6 (xi - X2) 2 J 

Again, as a becomes large x\ — > X2 and the effect of the defect on the singular solutions 
u and v disappears. However, the non- locality will require a large value of X\ — x 2 - 



11. Discussion 



The purpose of this investigation has been to discover to what extent the properties 
of jump-defects, originally described in certain relativistic integrable field theories, 
extend to non-relativistic systems of various kinds. In all cases, the defects are purely 
transmitting, in the sense that solitons passing through them may be delayed, converted 
to anti-solitons in some cases, or absorbed, but will never be reflected. In every case 
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examined, the integrable defect conditions investigated would constitute a Backlund 
transformation except that the spatial derivatives of the fields are frozen at the defect 
location. It is known that Backlund transformations for a given model are not unique 
and, particularly for the KdV and mKdV models, several different expressions are 
available in the literature. As Backlund transformations, all the available expressions 
are equivalent; but, as defect conditions, this does not appear always to be the case. For 
instance, for both the KdV and the mKdV equations local expressions, using u instead 
of p = u x , for a Backlund transformation are available [18]. However, it does not seem 
that all formulations of Backlund transformations may be used as defect conditions 
since it is not always possible to give a Lagrangian description for the ^-expressions for 
these models, and consequently to find suitable defect potentials. Moreover, to allow 
momentum conservation, which is a key feature of the jump-defect, the corresponding 
defect conditions for KdV and mKdV need to satisfy the additional relation listed at 
the end of section 9. Not all Backlund transformations seem to satisfy such a condition 
- for an example of one that does not, consider the Backlund transformation proposed 
in [21] for the KdV equation. 

It ought to be possible to modify appropriate Lax pairs for KdV and mKdV [24] 
to accommodate the jump-defects. However, because of the proliferation of higher 
derivatives in the expressions for the Backlund transformations this appears to be less 
than straightforward and discussions of these will be deferred. 

The nonlinear Schrodinger model seems to provide the nicest results. In this case, 
bound states have been found and it has been demonstrated using the two-soliton 
example that several solitons passing a defect will be affected by it independently 
of one another. If there are several jump-defects at different locations they too will 
act independently. It transpires that the strangest situation is illustrated by the KdV 
example. First, it does not seem possible to recover the single bulk soliton by taking a 
suitable limit of the defect parameter, and second the consequences of taking the limit 
\/3\ — > remains mysterious. In fact, the situation in which \/3\ < a, that is the defect 
parameter is less than the soliton speed, does not seem to be allowed since the incoming 
soliton at the left of the defect would become singular on the right of it. One resolution 
of this that has been suggested is that the soliton is trapped by the defect (as is clearly 
the case when \/3\ = a). 

It was noted in [13] that the jump-defects of the sine-Gordon model may themselves 
move and indeed scatter if there are several moving with different speeds. This aspect 
has not been considered in the present study though similar properties are expected, 
at least for NLS. The non-local nature of the defect conditions for the other models 
(KdV and mKdV) may pose a problem in this regard. It was also noted in [13] that the 
natural solution of the 'triangle relations' [14], 

s^(e)Xi(o 1 )Xi(o 2 ) = ^(e 2 )^(e 1 )s^ n (e), e = e 1 -e 2 , (n.i) 

where the roman labels are ±1, corresponding to soliton or anti-soliton, the greek labels 
are even integers, corresponding to the 'topological charge' carried by the defect, S^ n (<S)) 
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is the bulk scattering matrix for a pair of solitons, and c T^f a (6) is the transmission 
matrix for a soliton passing through an even-charged defect, is consistent with the 
Lagrangian description of the sine-Gordon jump-defect. In the present context, the 
NLS model appears to be the most suitable for a quantum investigation and also the 
most interesting due to the presence of bound states. It is well known that the quantum 
version of NLS is equivalent to a one dimensional multi-particle problem with 5-function 
pairwise interactions of equal strength. The problem in the bulk has been studied long 
ago, in the first place by means of the coordinate Bethe ansatz [27] (but see also [28] for 
a review and more references). It is therefore natural to try to incorporate the jump- 
defect within the N body quantum picture and to discover what manner of particle 
interaction is able to describe it. More recently, an N body problem of this kind with 
impurity (in a 'repulsive' regime where there would be no solitons in the corresponding 
NLS model), has been analysed in [29]. However, more investigations are necessary 
to see if there are any connections between that type of impurity and the jump-defect 
described in this article. 

One aspect of the story that remains frustrating is the absence of a physical 
model of an integrable jump-defect. It would be very interesting to discover a physical 
situation where the Backlund transformation plays a natural role. There are jump- 
defects commonly occurring naturally - modelling a dislocation in a material provides 
an example, as does modelling a shock front, or a fluid bore - where certain physical 
quantities are regarded as discontinuous, such as the fluid velocity on either side of the 
bore, yet others are continuous; and some conservation laws are preserved. However, 
there does not yet appear to be a nonlinear, and integrable, example of the kind we 
are discussing. If there were, it might offer the possibility of controlling solitons, which 
might then be put to use (see [30] for a suggestion). 
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